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towards on-demand generation of photon-number states 



I. Dotsenko 1 - 2 ,* M. Mirrahimi 3 , M. Brune 1 , S. Haroche 1 ' 2 , J.-M. Raimond 1 , and P. Rouchon 4 

Laboratoire Kastler Brossel, Ecole Normale Superieure, CNRS, 
Universite P. et M. Curie, 24 rue Lhomond, F-75231 Paris Cedex 05, France 
2 College de France, 11 Place Marcelin Berthelot, F-75231 Paris Cedex 05, France 
3 INRIA Rocquencourt, Domaine de Vouceau, B.P. 105, 78153 Le Chesnay Cedex, France 
4 Mines ParisTech, Centre Automatique et Systemes, Mathematiques et Systemes, 
60 Bd. Saint-Michel, 75272 Paris Cedex 06, France 
(Dated: May 1, 2009) 

We propose a quantum feedback scheme for the preparation and protection of photon number 
states of light trapped in a high-Q microwave cavity. A quantum non-demolition measurement 
of the cavity field provides information on the photon number distribution. The feedback loop is 
closed by injecting into the cavity a coherent pulse adjusted to increase the probability of the target 
photon number. The efficiency and reliability of the closed-loop state stabilization is assessed by 
quantum Monte-Carlo simulations. We show that, in realistic experimental conditions, Fock states 
are efficiently produced and protected against decoherence. 

PACS numbers: 42.50.Dv, 02.30.Yy, 42.50.Pq 



I. INTRODUCTION 

It is now possible to realize ideal quantum measure- 
ments on individual quantum objects, for instance atoms, 
ions and photons. Beyond being a tool to monitor the 
system's evolution, projective quantum measurements 
can be used to prepare it in specific quantum states. For 
instance, an ideal quantum non-demolition (QND) mea- 
surement of a field's photon number projects it onto a 
photon-number (Fock) state [1]. However, due to the 
basic quantum indetermination of the measurement out- 
come, measurement-induced state generation is not de- 
terministic. Quantum feedback control techniques [2] 
make it possible to overcome this limitation and to pro- 
duce quantum states on demand. These techniques gen- 
erally combine weak quantum measurements with a real- 
time correction of the system's state depending on the 
classical information extracted from the measurements. 
Beyond preparation of specific states, these feedback 
schemes can also protect them from decoherence, result- 
ing from the coupling of the system with its environment. 

In this paper, we propose a quantum feedback scheme 
for the on-demand preparation of Fock states stored in 
a high-quality superconducting microwave cavity and for 
their protection against decoherence. This scheme is de- 
signed to operate with an existing cavity-QED set-up [3]. 
The feedback loop uses three steps. We first extract infor- 
mation on the photon number distribution with a single 
circular Rydberg atom by a QND process. This atom 
interacts dispersively with the cavity field and does not 
exchange energy with it, but experiences a light shift pro- 
portional to the photon number. The measurement of 
this shift, with the help of a Ramsey interferometer, pro- 
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vides information on the field intensity. It modifies the 
field state accordingly through the quantum projection. 
This information is, in the second step, used to estimate 
the new cavity field state through a quantum filtering 
process [4, 5]. In the third step, we correct the field 
state, using a coherent field pulse injected in the cavity. 
The feedback law used to calculate the amplitude of this 
pulse is chosen by adapting to this discrete situation the 
Lyapunov-based techniques proposed in Ref. [6]. These 
three steps represent the three basic components of any 
feedback loop: a sensor, a controller and an actuator. 
In contrast to classical closed-loop systems, we use here 
a quantum sensor. However, both controller and actua- 
tor remain classical. By iterating the feedback loop, we 
steer the cavity towards any target Fock state. We also 
efficiently stabilize it against cavity decay. 

Recently, a similar feedback scheme to generate Fock 
states of light in an optical cavity has been proposed [7] . 
It is based on a continuous monitoring of the mean num- 
ber of photons. Instead, we propose here to use a discrete 
QND measurement followed by a quantum filtering pro- 
cess providing complete information on the field's den- 
sity matrix. This precise knowledge of the field's state 
gives us a much better insight into the feedback action. 
It could also be an asset in other feedback schemes re- 
quiring, for instance, information on the field's phase in 
addition to its intensity. 

The paper is organized as follows. In Sec. II, we de- 
scribe the elements of the set-up and give realistic experi- 
mental parameters. We present the quantum-mechanical 
operators describing the evolution of the cavity state un- 
der measurement, decoherence and pulse injection. Sec- 
tion III is devoted to a detailed analysis of the quantum 
filter and of the feedback law in an ideal situation. We 
first present the main elements of the feedback loop and 
describe the tuning of the controller gain. We give qual- 
itative arguments showing that the proposed strategy is 
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Figure 1: Proposed quantum feedback scheme adapted to a 
microwave cavity QED set-up. C: high-Q microwave cavity, 
B: box producing Rydberg atoms, R\ and ife: low-Q Ram- 
sey cavities, D: atomic field-ionization detector, S and S': 
microwave sources coupled to C and R's, respectively. In a 
quantum filtering process, a real-time control system analyzes 
the results of QND measurements of the cavity field and com- 
putes the amplitude of a control injection pulse. 



stable. The detailed mathematical proof of convergence 
and stability, relying on stochastic Lyapunov techniques 
will be published elsewhere [8]. We finally present quan- 
tum Monte Carlo simulations of closed-loop trajectories 
of the cavity field exhibiting the feedback performances. 
In section IV, we take into account the known experi- 
mental imperfections of the existing set-up. We modify 
accordingly the feedback algorithm and present extensive 
simulations of its operation. We conclude in Sec. V. 



II. EXPERIMENTAL SET-UP 

Our quantum feedback algorithm is designed for the 
ENS microwave cavity QED set-up [3, 9]. Its compo- 
nents are depicted in Fig. 1. The microwave field to be 
controlled is confined in an ultra-high Q superconducting 
cavity C (damping time T cav = 0.13 s). Rydberg atoms, 
flying one by one at 250 m/s across the cavity mode and 
dispersively interacting with its field, are used as QND 
probes of light [1]. They are prepared in a superposi- 
tion of two circular Rydberg states |e) and \g) (principal 
quantum numbers 51 and 50, respectively) in the low-Q 
cavity Ri, driven by the source S' , in which they expe- 
rience a resonant 7r/2-pulsc. In the Bloch sphere repre- 
sentation, the spin corresponding to the two-level atomic 
system then points along the Ox direction in the equa- 
torial plane (states |e) and \g) correspond to north and 
south poles of the sphere, respectively; we use a frame 
rotating at the atomic frequency). 

The light shifts experienced by the non-resonant \g) — > 
|e) atomic transition in the cavity field result in a phase 
shift for the atomic state superposition. At the cavity 
exit, the atomic spin points along a direction in the equa- 



torial plane at an angle $(n) with the Ox axis, correlated 
to the photon number n. The dephasing angle $(n) can 
be controlled by adjusting the atom-cavity detuning S 
and the interaction time. In the large atom-cavity detun- 
ing regime, it is a linear function of n: $(n) ~ (Sl^/4S)n, 
where rio/2ir = 49 kHz is the vacuum Rabi frequency. 
For the intermediate detuning values used in the experi- 
ments, it is a more complex, albeit perfectly known grow- 
ing function of n. A second ir/2 Ramsey pulse in R 2 
(phase </>r with respect to that of the pulse in Ri) fol- 
lowed by the atomic detection (in the {|e) , \g)} basis) 
by the detector D amounts to a detection of the atomic 
spin along an axis at an angle c/)r, with Ox. It provides 
information on the cavity field intensity. 

Atoms are sent in the set-up at typical 250 p,s time in- 
terval, much shorter than the cavity damping time. Note 
that such a macroscopic time interval is well adapted to 
elaborate feedback strategies since we have ample time 
to compute the state estimator and the feedback law be- 
tween two atomic detections. When no feedback action 
is performed, the information provided by a few tens of 
atoms results in a measurement of the dephasing angle 
$(n) and, hence, in a projective QND measurement of 
the photon number n [1, 10]. 

We are interested here instead in the ambiguous infor- 
mation provided by a single atomic detection. Detection 
of the atomic state in D projects the field, described ini- 
tially by the density matrix p, onto a new state p P roj- 
Depending on the detected atomic state, \s) = |e) or \g), 
the back-action of the quantum measurement on the field 
is described by 

M sP Ml 

Pproj — / ,\ ) i 1 ) 

Tr[M s pMj) 
where the operators M s are given by 

Ma=cos (*L±im), (2a) 

Mt= , in (*siSl). (2b) 

Here, N = a^a is the photon number operator with a 
and or the photon annihilation and creation operators, 
respectively. The measurement operators M s are diag- 
onal in a photon- number basis (as is N). They thus 
preserve Fock states, illustrating the QND nature of the 
measurement. The projected state (1) is normalized by 
the probability P s = Tr (M s pMJ) of detecting the atom 
in state \s). 

The cavity field can also be manipulated by injecting 
into the mode a coherent field pulse generated by the 
resonant microwave source S. Its action is described by 
the displacement operator D{a) = exp(aa^ — a* a), where 
a is the complex amplitude of the injected field. The 
cavity field after displacement thus reads 

Pdisp = D(a) pD{-a). (3) 



3 



A proper analysis of the feedback scheme requires to 
take into account all known imperfections of the experi- 
mental set-up. Besides projective measurement (1) and 
coherent evolution (3), the cavity state also evolves due 
to decoherence. The field is coupled to a reservoir at 
non-zero temperature (T = 0.8 K), and its dynamics is 
described by the master equation [11] 

^ = Lp = - - (1 + n t h)(a t ap + paa* - 2apa f ) 

— 2 n th(o,a ! p + pa 1 a — 2a) pa), (4) 

where k = 1/T cav is the cavity decay rate and nth = 0.05 
the equilibrium thermal photon number. 

Moreover, the circular Rydberg state preparation is a 
non-deterministic Poisson process. We perform a pulsed 
excitation from a continuous thermal beam of Rubidium 
atoms, preparing atomic samples at a T a = 85 /is time 
interval. The mean number of circular atoms per sample 
is 7? a sa 0.3 (leading to an average time between atoms 
of 250 ps). The atomic time of flight between C and 
D is 350 (Lts, corresponding to a delay of d = 4 atomic 
samples. The probability for detecting an atom present 
in a sample is rjd ~ 0.8 due to finite detection efficiency 
in D. Finally, non-ideal atomic state resolution of the 
ficld-ionization detector and finite contrast of the Ramsey 
interferometer introduce a probability of erroneous state 
assignation of i]f m 0.1. 



III. IDEALIZED EXPERIMENT 

In order to describe the essential elements of the feed- 
back scheme in a simple context, we first consider in 
this Section an idealized experiment with no cavity decay 
(k = 0), a deterministic atomic preparation (?y a = 1), no 
detection delay (d = 0), and a perfect atomic detection 
(rid = 1 and r?f = 0). 



A. Quantum filter 

The quantum filtering procedure [4, 5] provides us with 
an estimate of the field state by using all available infor- 
mation. Right after detection of atom number k and 
before injection of the corresponding control, it includes 
the initial state of the field po , information obtained from 
all k atoms detected so far (labelled i = 1 . . . k) and all 
coherent field injections performed in the former feed- 
back loops (1 to k — 1). Each of the elementary pro- 
cesses (atomic detections and displacements) can be rep- 
resented by super-operators acting on the field's density 
matrix. Let us note M, that associated to the detection 
of atom i and D i that corresponding to the displace- 
ment performed in the i th iteration of the feedback loop. 



Therefore, after detecting atom k, the state of the field is 

Pk = M k (f[-D i Mi\ Po- (5) 

Here and in the following, all super-operator products are 
ordered by decreasing indices. Note that the expression 
(5) can be computed iteratively from the recurrence 

p k +i = M k+1 T> k p k , (6) 

setting D = 1. By applying this filter recursively at 
each cycle, we get a reliable real-time estimate of the 
field state. 

Based on Eqs. (1) and (3), the super-operators are de- 
fined as 

M SiP Mj 

M iP = 7 TV' ( ? ) 

Tr [M Si pMlj 
D i/9 = D(ai)pD(-ai). (8) 

Here, Sj = g or e depending on the outcome of the i th 
atomic measurement. The control amplitude at is ad- 
justed after each atom detection according to the feed- 
back law described in the next subsection. 

We choose to prepare initially the cavity in a coherent 
state with a mean photon number equal to the target 
value n tag . The phase of this field is used as a refer- 
ence. As a consequence, its amplitude ao = y/n tas is 
real. The phase 4>r an d the operator $(AT), defining the 
measurement super-operators M, through (2), should be 
adjusted to optimize the quantum filter (6). On the one 
hand, the mean dephasing angle per photon </> (estimated 
in the useful range of n values) should be large in order to 
resolve adjacent photon numbers. On the other hand, the 
measurement should lift the ambiguity between a large 
range of possible photon numbers from to some n max , 
determined by the spread of the initial coherent field. 
Due to the oscillating nature of the measurement opera- 
tors Ms, this sets a maximum value of <p of the order of 
7r/n max . 

Finally, we need to optimally distinguish n tag from 
n tag ± 1. This is achieved by setting the Ramsey phase 
0r so that (</>r + $(n taB ))/2 = ir/4. This setting cor- 
responds to a maximal variation of the atomic state de- 
tection probabilities around n tag . In other words, we set 
the atomic Ramsey interferometer at mid-fringe when the 
cavity contains n tag photons. 



B. Feedback control 

The distance between p k and the target Fock state de- 
scribed by a density operator p tag = |n tag ) (n tag | can be 
conveniently defined as 

V(p k ) = l-Tr(p k p t J = l-F(p k ), (9) 
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where F(p) = Tr(pp tag ) is the fidelity of a state p with 
respect to the target state p tag . The function V(p) rep- 
resents a natural choice for the Lyapunov function [6] 
used to study the stability properties of the quantum con- 
troller. This study is presented in more detail in Ref. [8]. 

Let us denote the cavity state after the control dis- 
placement in the fc th iteration of the feedback loop as 
p k = D k p k . The optimal convergence of the feedback 
procedure is obtained by finding at each cycle the injec- 
tion amplitude a k which maximizes the fidelity F{p k ). 
This search can be performed in a straightforward way by 
numerical optimization. However, this approach is time- 
consuming and could not be realistically implemented 
with existing real-time data analysis systems in the few 
tens of microseconds time interval between two atomic 
detections. We overcome this difficulty by developing 
a faster, analytical calculation of a locally optimal dis- 
placement a k . Since the measurement operators, the ini- 
tial density matrix and the projector on the target Fock 
state are all real operators, we can reasonably choose a k 
to be real. The field density matrix then remains always 
real. 

In the limit of small a, the Baker-Campbell-Hausdorff 
formula [12] yields the following approximation for a dis- 
placed state (3): 

p = D(a)pD(—a) «p - a[p, — a] 

2 

+ y[by -a], a f -a}. (10) 
Using (8), we therefore get for small a k 
F{p k ) = F(p k ) - a k Tr ([p k , aJ — a]p t3 , g ) 

+ &Tr([\p k ,a!-a],<f-a]p u .) . (11) 
Choosing the feedback amplitude as 

a k = -CiTr ([p k ,al -a]p tas ) , (12) 
with a gain c\ > small enough ensures that 

F(p k ) > F(p k ) (13) 

as soon as the trace in (12) is strictly positive. Fur- 
thermore, since [p tag ,M g ] = [p tag ,M e ] = and M^M g + 
M\M e = 1, the conditional expectation of F(p k+ i) 
knowing p k is given by 

p t ( P^ M a~Pk M t \ 
^ (F{Pk+i) | Pk) = P g ,kTr I p— I 

+ P e ,^( A ^ Ml )=F(p fe ). (14) 

Therefore 

E(F(p k+1 ) | p k ) >E(F( Pk+1 ) | p k ) = F(f>k) (15) 



and, as a result, the expectation value of F(p k ) is a non- 
decreasing function of k: 

l>E(F(p k+1 ))>E{F(p k ))>0. (16) 

We have only shown so far that F(p k ) is increasing on 
the average. It does not mean that, for all individual re- 
alizations of the loop sequence, F(p k ) converges to 1, its 
maximum value reached when p is equal to p tag . Indeed, 
the feedback law (12) does not prevent the convergence 
towards other Fock states, since a k = whenever p k is 
the projector on any photon number state. A careful 
analysis shows that, in each realization, F(p k ) converges 
cither to 1 or, with some small probability, to zero. To 
overcome this spurious attraction towards wrong photon 
number states, we propose, as in [6], to modify the feed- 
back law (12) by applying a constant injection "kick" as 
soon as the cavity state significantly deviates from the 
target: 

JciTr ([p tag ,al-a}pk) if F{p k ) > e, 
(17) 
c 2 sign (n tag - n k ) if F(p k ) < e. 

Here, n k is the mean photon number in the current state 
Pk, c-2 > is a constant kick amplitude. The kick condi- 
tion is set by e with 1 > 0. In this way, by breaking 
attraction towards other Fock states, we make sure that 
the cavity state converges towards the target one. 

The controller gain c\ in (17) must be tuned to maxi- 
mize fidelity F(p k ) at each sampling time k. Up to third- 
order terms in (p k — p tag ), Eq. (11) yields 



F(p k ) = F(p k ) + 




By maximizing this expression, we find that the gain 

ci =Tr([ (0tag , a t_ a ]2)- 1 = (4n tag +2)- 1 (19) 

provides the fastest convergence speed in the vicinity of 
p tag . A detailed mathematical analysis and a convergence 
proof of this feedback-scheme will be given in [8]. 

C. Simulations 

To assess the performance of the proposed feedback 
scheme, we have performed quantum Monte-Carlo nu- 
merical simulations. Figure 2(a)-(c) presents a single 
quantum closed-loop trajectory in the ideal measurement 
setting, with a target Fock state |n t a g ) = |3). The Hilbcrt 
space is limited to n max = 9 photons. The atom-cavity de- 
tuning is set to 5/2ir = 238 kHz. The mean dephasing an- 
gle per photon around n tag is then about tt/7. Intuitively, 
the Ramsey phase </>r should be set at the mid-fringe 
setting for n tag photons, i.e., (0 R)O + $(n tag ))/2 = 7r/4. 
We have numerically observed that the convergence of 
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Figure 2: A single closed-loop quantum trajectory of an ideal 
measurement with |n tag ) = |3). (a) Atomic states detected 
at each feedback cycle, (b) Control injections at- (c) Evo- 
lution of photon number probabilities P(n < n ta& ), P(n tag ) 
and P(n > n tag ) shown as dashed-dot blue, solid green and 
dashed red curve, respectively, (d) Average over 10 4 closed- 
loop quantum trajectories. 



the process is faster when the Ramsey phase alternates 
between four values in successive feedback cycles: </>r,o, 
(0R,o + tr), 0r,o and (<Ar,o — cr)- The phase excursion 
or is set at 0.69 rad. The kick amplitude ci is 0.1 with 
a kick zone defined by e = 0.1. The initial cavity field 
po is the coherent state with 3 photons on the average, 
Po = D (V3) |0) (0\D(-VS). 

The upper trace in Fig. 2 shows a sequence of detected 
atomic state. The corresponding control inputs are 
plotted in the next panel. Large at the beginning, cv^ 
rapidly converges towards zero. Figure 2(c) shows the 
probability for the number of photons in the cavity field 
to be smaller, equal or larger than n tas (dashed-dot blue, 
solid green and dashed red curve, respectively). For this 
particular quantum trajectory, we observe that p con- 
verges to p tag in less than 30 feedback cycles. 

Figure 2(d) presents an ensemble average over 10 4 real- 
izations of the same numerical feedback experiment. The 
fidelity F{pk) reaches values above 0.99 after 140 cycles. 
After 20 cycles only, about 80% of the trajectories have 
converged to the target state, illustrating the efficiency 
of the feedback method. 



IV. REALISTIC EXPERIMENT 

So far, we have neglected experimental imperfections. 
In this Section, we take into account all known imperfec- 
tions of the present experimental set-up: finite cavity life- 
time, Poisson distribution of the atom number in atomic 
samples, non-ideal efficiency and state-selectivity of the 
detector and, finally, the finite delay between atom-cavity 
interaction and atomic detection. Below, we introduce 
modifications of the quantum filter (6), which make the 
process efficient in spite of these imperfections. 



A. Modified quantum filter 

The quantum filtering process must provide us with a 
density matrix describing all our knowledge of the cavity 
field at a given time. We chose here this time to be right 
before the injection of the actuator pulse in the k th iter- 
ation of the loop, after detection of the k th atom. At this 
time, our knowledge includes all atomic measurements 
so far (from 1 to fc), all displacements performed and 
the known relaxation of the cavity field in between these 
events. It must also include the influence of the d atomic 
samples which have interacted with the cavity, but are 
still flying towards the detector when we record atomic 
sample k. All this information can, as above, be rep- 
resented by super-operators acting on the field's density 
matrix. 

Let us note P s the super-operator describing the de- 
tection of an atomic sample. Each detection has three 
possible outcomes, labelled s = e, g or u (atom detected 
in e, g or no detected atom at all). The super-operator 
N describes the action of one of the d atomic samples 
flying from C to D. The displacement super-operator 
D is given by (3) and the field relaxation during the 
time interval between atomic samples is described by T . 
The exact expressions of all these super-operators will be 
given later. 

Therefore, the state of the cavity field before injec- 
tion is 



Pk — F 2 ,k Fi.fc Po, 



where 



fc-i 



F lfc = TP Sfc ]JU i_ d TP 3i 



(20) 



(21) 



includes the information gathered by the first k sample 
detections and 



k+d-l 
i—k 



(22) 



includes the influence of the d in-flight samples. Note 
that due to the atomic propagation delay between C and 
D, the displacement amplitude injected in C right after 
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atomic sample i has interacted with it is czi-d-, computed 
in the (i — cQ-th iteration of the loop. Hence, the first d 
atomic samples cross the cavity before any displacement. 
This is formally taken into account by setting D i to unity 
for non-positive indices. If d = 0, i.e. for no delay in the 
detection process, the empty product in (22) equals by 
convention the identity operator. Consequently, in the 
case of the ideal experimental parameters, the expression 
(20) naturally reduces to (5). 

The density operator p^ is used in the feedback law 
(17) to calculate the control injection ctk- The product 
Fi : fc can be easily computed recursively, as in the ideal 
case. However, we must recalculate F-2,k at each feedback 
cycle. 

We now give the explicit expression of all super- 
operators considered so far. The measurement process 
P s takes into account detector's imperfections and non- 
ideal atomic source. If the detector clicks and indicates 
state \s) = \g) or |e), the new quantum state is a sta- 
tistical weighted mixture of two projected states. The 
strongest weight corresponds to the projection according 
to the recorded outcome. With a smaller probability, the 
atomic detection has been erroneous. The atom was, at 
detection time, in the state \s), opposite to \s), and the 
field has been projected accordingly. Thus, due to the 
imperfections of state discrimination alone, we get 

P (s=e / 9 ) P = (1 - Pt)M sP + PMsP ■ (23) 

The weights of the two projected states are given 
by the conditional probability of a wrong detection 
knowing that the detector has clicked in |s): Pf = 
m P-/ ((1 - m )P s + m P-) with P s = Tr (M sP Mt). 

If the detector does not click (outcome s = u), either 
there was no atom in the sample or an atom was present 
but has not been detected. Therefore, the estimated state 
is a mixture of the unperturbed cavity field and of the two 
projected states corresponding to the two possible states 
of the undetected atom (we assume here that these two 
states are equiprobable) . We thus get 

P(«=u) P = {^-Pu)p + Pu (M g pMl+M e pMl) . (24) 

The weights in this mixture are set by the conditional 
probability to have an undetected atom in a sample: 
Pu = ?7a(l - rjd)/ (1 - IhVd)- 

The interaction with a not-yet-detected sample is de- 
scribed by the super-operator N. It is given by (24) with 
the conditional probability P u replaced by the atomic 
sample occupation probability 7y a , since the measurement 
has not been performed yet, resulting in 

Np = (l-r/ a )p + r h (M gP Ml+M eP Mt) . (25) 

The super-operator T describes the evolution due to 
the cavity field relaxation during the time interval T a 
between two atomic samples. Using (4) in the approxi- 
mation of small time interval, T" a ^;T cav , we get 

Tp=(l + T a L)p. (26) 



B. Simulations 

We have performed extensive simulations, similar to 
those of the idealized measurement in Fig. 2, with the 
feedback algorithm adapted to the realistic experimental 
parameters, given in Sec. II. We use Monte-Carlo sim- 
ulations to follow the state of the cavity along individ- 
ual realizations of the experiment. The field relaxation 
is taken into account by simulating random quantum 
jumps, whose average effect is described by Eq. (26) [13]. 
For each atomic sample, we choose randomly the number 
of atoms (0 or 1) and the real atomic state (|e) or \g)) 
which would be detected by the ideal detector. These re- 
sults determine uniquely the subsequent field state p™ 1 . 
We then take into account the detection imperfections 
to simulate the outcomes of the realistic measurement, 
which arc used in the feedback loop for computing pj~. 
Note that /o real is not accessible in the experiment. Here, 
we use it to simulate the evolution of the field and, by 
comparing it with the state estimated in the quantum fil- 
tering process, to characterize the feedback performance. 
In the ideal case, considered in Sec. Ill, the estimated 
state pk always coincide with the real state of the field 
PT l - 

Figure 3(a)-(d) shows a typical quantum trajectory, il- 
lustrating the robustness of our feedback scheme at the 
present level of imperfections. The first difference with 
the ideal case is the presence of many feedback cycles 
with no detected atoms. Nevertheless, after about k — 
100 feedback cycles (each of 85 ps duration), correspond- 
ing to approximatively krj a ss 30 atoms and kr/^rj^ f=a 24 
detector clicks, the cavity state successfully converges to 
the target Fock state |n taB ) = |3). 

The other striking feature is the presence of sudden 
photon jumps due to the limited cavity lifetime. In this 
particular trajectory, a photon loss occurred at feedback 
cycle number k « 430 (at about 37 ms). One sees that 
about 60 cycles (about 5 ms) are required to detect the 
jump before applying a large control displacement and 
about 60 more cycles to restore the target photon number 
in the cavity field. 

Figure 3(e) shows the same type of ensemble average as 
presented in Fig. 2(d). We observe a convergence towards 
a fidelity of 0.63 after several hundred cycles resulting in 
about 100 detector clicks. This convergence, given as 
the number of actually detected atoms, is similar to that 
observed in the ideal case. The main difference between 
the two traces is the reduced asymptotic value of the 
fidelity of the prepared state: F = 0.63 instead of F = 1 
in the ideal case. This reduction is mainly due to the 
cavity decay, which results in uncontrolled jumps of the 
photon number. The observed fidelity indicates that for a 
typical closed-loop trajectory the field stays in the state 
1 3) about 63% of time. Since the lifetime of |n tag ) is 
given by T cav /((1 + n t h)n ta& + n th (n + 1)) = 45 ms [10], 
we estimate that it takes on the average 26 ms to restore 
the initial Fock state after a sudden jump. 

In order to observe in more detail the recovery of the 
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Figure 3: A single closed-loop quantum trajectory for realistic 
experimental parameters with |n tag ) = |3). (a) Sequence of 
detected atomic states, (b) Control injections a^. (c) Evolu- 
tion of the photon number probabilities P(n<n tag ), P(n tag ) 
and P(n>n tag ) (dashed-dot blue, solid green and dashed red 
curve, respectively) of the field state p" ai deduced from the 
Monte-Carlo simulation. The thin black line gives the esti- 
mated state fidelity F(pk). (d) Photon number probabilities 
for n = to 9 in p real . For better clarity, results of only one 
out of each 10 feedback cycles are shown, (e) Average over 
10 4 closed-loop quantum trajectories. If the number of quan- 
tum trajectories is large enough (as is the case here), the av- 
eraged photon number probabilities calculated from the real 
field states p rcal coincide with those averaged over the esti- 
mated states. 




10 20 30 
Time [ms] 

Figure 4: Feedback performance in the presence of a sud- 
den photon number loss. The fidelity of the actual cavity 
state p™ 1 (blue solid curve) and of the state estimated by the 
quantum filtering process pk (red dashed curve) are defined 
with respect to the target Fock state and are obtained by av- 
eraging over 10 4 quantum trajectories. In each trajectory, the 
time origin is shifted to the time of a quantum jump. 



system after a sudden photon jump, wc recorded 10 4 
quantum trajectories showing sudden photon losses at 
different times. Before averaging the trajectories, we shift 
their individual time origins to the time of a jump. Fig- 
ure 4 shows the evolution of the fidelity of the actual 
(solid curve) and estimated (dashed curve) cavity states 
with respect to p tag . As calculated before, it takes about 
20 — 30 ms to restore the initial photon number state 
after a sudden photon loss. 

The maximum fidelity of the quantum feedback, seen 
in Figs. 3-4 and corresponding to the ensemble of many 
realizations, shows how well a chosen photon number 
state can be preserved on average from its unavoidable 
decay. However, since in each feedback cycle we acquire 
actual information on the field, the fidelity of the deter- 
ministic generation of Fock states in a single realization 
can be much higher. To reliably produce the target state 
p tag we keep the feedback process running until we de- 
tect its successful convergence to p tag with an estimated 
fidelity better than F com . Figure 5(a) shows a histogram 
of convergence times for F conv = 95 % obtained from 10 4 
quantum trajectories. The cumulative distribution, pre- 
sented by a solid line, shows that after 20 ms the proba- 
bility of a successful feedback outcome reaches 50 % and 
that it exceeds 90 % after 85 ms of the feedback opera- 
tion. The inset in Fig. 5(a) shows the average density 
matrix of converged states with the population in |n tag ) 
of about 95%. 

Our feedback scheme can be tuned to generate any 
small photon number states. Figure 5(b) illustrates the 
convergence for several Fock states from |n tag ) = |1) to 
|5). By increasing n tag , the initial photon number distri- 
bution gets broader as the photon number lifetime gets 
shorter. Therefore, the feedback performance slowly de- 
grades, resulting in the increased convergence time for 
the same required fidelity F conv . 
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Figure 5: Feedback convergence, (a) Histogram of conver- 
gence times with F conv = 0.95. The solid line gives the 
probability of the feedback convergence after a given time. 
The inset shows the average density matrix of the converged 
states, (b) Convergence for several different Fock states. 
These curves result from an average over 10 4 quantum trajec- 
tories. 



C. Real-time computation 

Implementation of this quantum feedback scheme into 
a real experiment requires fast real-time analysis of the 
measured data, faster than duration T a of one feedback 
cycle. We have tested a real-time processing system 
ADwin (Jager Messtechnik) allowing for fast data ac- 
quisition and analysis in combination with complex ex- 
periment control. After careful optimization of all feed- 
back computations, the execution of one typical feed- 
back cycle of the proposed scheme requires about 4000 
floating-point operations. Our preliminary tests have re- 
vealed that the ADwin system needs for this task about 



30 fis < T a = 85 //s. We are now working on the integra- 
tion of this control system into our cavity QED set-up. 



V. CONCLUSION 

We have presented a quantum feedback protocol de- 
signed to deterministically generate small photon number 
states of a trapped microwave field. The reliable estima- 
tion of the cavity state at each feedback cycle allows us 
to follow in real time the quantum jumps of the photon 
number and then to efficiently compensate them, thus 
protecting Fock states against decoherence. 

The two main components of the feedback algorithm 
are the quantum filter, which estimates the actual state 
of the field based on the outcomes of QND measurements, 
and the feedback law computing the amplitude of a cor- 
rection microwave pulse whose injection into the cavity 
mode maximizes the field fidelity with respect to a de- 
sired Fock state. Quantum Monte-Carlo simulations of 
the QND measurements and the quantum feedback re- 
sponse demonstrate the high reliability of our closed-loop 
scheme even in the presence of realistic experimental im- 
perfections. This scheme is also robust with respect to 
an imprecise knowledge of the initial state, since the re- 
peated measurements of the state corrects initial wrong 
or incomplete information. Convergence and stability 
proofs of these feedback schemes will be given elsewhere. 
For a first set of mathematical results, see Ref. [8]. 

The presented scheme can be extended to include 
the Ramsey phase </>r and the spin dephasing $(iV) 
as additional control parameters of the feedback loop. 
This should help to optimize the QND measurement 
according to our estimate of the actual field state and 
thus to acquire more information on its dynamics. These 
studies arc under the way. 
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